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ABSTRACT 

Aims. We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object 
(YSO), at very high spatial resolution (10-100 AU). 

Methods. We performed milli-arcsecond polarimetric observations of the strong CH 3 OH maser emission observed in the vicinity of 
an 0-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH 3 OH masing gas previously 
measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction 
of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. 
Results. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a 
massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, 
where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2° AU“'). Overall, the velocity field 
vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30°. We interpret this 
finding as the contribution of a turbulent velocity field of about 3.5 km s“', responsible for braking up the alignment between the 
velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk 
plane, with an average speed of 7 km s“'. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance 
of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 
500-600 AU, and inward for shorter radii. 
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1. Introduction 


The role of magnetic fields in regulating the gas dynamics in the 
vicinity of growing, massive, young stellar objects (YSOs) is still 
a matter of debate (e.g., [Crutcher et al.|2010HZhang et d? 2014| l. 
Recent magneto-hydrodynamics (MHD) models, simulating the 
build-up of massive protostars in the inner few lOOOAU, have 
shown that magnetic fields may contribute significantly (1) to 
the degree of outflow collimation, and (2) to stabilizing both 
Keplerian and sub-Keplerian disks against fragmentation (e.g., 
Seifried et al.|20ir |2012| l. In this context. Very Long Baseline 
Interferometry (VLBI) observations of maser emission, arising 
within a few 1000 AU from massive YSOs, allow us to determine 
both the velocity distribution and the magnetic field configura¬ 
tion close to the accreting protostar (e.g.,|Sanna et al.|2010a H 


[Goddi et al.||201l] [Moscadelli et al.||201 lUSurcis et al.||2015| l. 

This gives us the unique chance to investigate, at a 10-100 AU 
scale, whether or not the magnetic field influences the gas kine¬ 
matics. 

G023.01 -00.41 is a lumin ous star-forming region of about 
4 X 10"^ Lq ( |Sanna et al.||2014 l, located at a trigonometric dis¬ 
tance of 4.6 kpc ( [Brunthaler et al.||2009| l. This star forming site 
harbors a flattened, hot molecular core (HMC) which is centered 


on an active site of strong maser and radio continuum emis¬ 
sion ( jSanna et al.|2010b their Fig. 4). The kinematics of warm 
(200 K) gas in the inner 3000 AU, as traced with CH3CN and 
thermal CH3OH lines, shows the composition of two, orthogo¬ 
nal, velocity fields ( jSanna et al.|2014 their Fig. 3). The velocity 
component which dominates at larger scales is aligned with the 
axis of a collimated bipolar outflow, traced progressively away 
from the HMC center with SiO and CO gas emission. Since the 
outflow emission is almost perpendicular to the line of sight, any 
associated disk should be seen edge-on, which makes this object 
an excellent target to study the gas dynamics in the vicinity of 
an O-type YSO. Furthermore, the 3D gas kinematics revealed 


by the CH3OH masers shows a funnel-like morphology (Sanna 
et al.|2010b their Fig. 6), which was best interpreted as the base 


of the outflow cavity (or the surface of a flared disk) with a size 
between 1000 and 2000 AU. 


With this in mind, we decided to use the synergy between 
maser proper motions and polarization measurements, targeting 
the rich CH3OH maser spectrum observed in G023.01-00.41, 
to investigate whether magnetic fields may be actively driving 
the circumstellar gas motion around a massive YSO. That can 
be assessed by quantifying whether a correlation exists between 
the orientation of the velocity and polarization vectors locally, as 
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measured for individual CH 3 OH masing cloudlets on scales of a 
few AU. In order to compare the magnetic field orientation with 


the velocity held previously measured by Sanna et al. ( | 2010 b| l, 
we conducted polarimetric observations of the 6.7 GHz CH 3 OH 
masers towards G023.01-00.41 with the European VLBI Net¬ 
work (evnQ. 


2. Observations and Calibration 

We employed the EVN to observe in full polarization mode the 
5i - 6 oA^ CH 3 OH maser transition, at the rest frequency of 
6668.519 MHz, toward G023.01-00.41. The observations were 
conducted under program ES067 on 2011 May 29. We made 
use of a single frequency setup to obtain both a high spectral 
sampling (0.98 kHz) of the maser lines, and a bandwidth large 
enough (2 MHz) to accurately measure the continuum emission 
of the calibrator, J 2202-1-4216. This calibrator served both as a 
fringe hnder and polarization calibrator, and was observed ev¬ 
ery 45 min to properly calibrate the polarization leakage. Since 
many maser features are expected to be linearly polarized at a 
level of about 1 %, to reach a conservative detection above 5cr 
over half of the maser cloudlets previously detected (with peak 
intensities >3 Jybeam^^), we spent about 4.5 hours on-source. 
The EVN data were processed with the SEXC software corre¬ 
lator ( Keimpema et al.|[2015| l at the Joint Institute for VLBI in 
Europe by using an averaging time of 2 s. The single-dish spec¬ 
trum of the CH 3 OH maser emission toward G023.01-00.41 is 
plotted in Eig .[2 

Data were reduced with the NRAO Astronomical Image Pro¬ 
cessing System (AIPS). We mapped the CH 3 OH maser distri¬ 
bution with a (robust 0) beam size of 11 mas x 4 mas, achiev¬ 
ing a thermal noise of 5mJybeam“^ To calibrate the system¬ 
atic rotation of the linear polarization angle Ct^poi) in the EVN 
dataset, we compared the EVN measurement of^poi obtained on 
J 2202-1-4216 with two, consecutive, VLA polarimetric observa¬ 
tions of the same calibratoij^ bracketing our VLBI observations 
(on 2011 April 30, and 2012 Eebruary 3). The linear polarization 
angle of J 2202-1-4216 estimated with the VLA remained nearly 
constant with an average value of-31° + l° (position angles, e.g., 
;t'poi, are measured east of north, unless otherwise stated). There¬ 
fore, the Xpo\ measurements obtained with the EVN dataset are 
affected by a systematic uncertainty of no more than a few de¬ 
grees. The uncertainty of;t'poi due to thermal noise was obtained 
from the relative error of the polarization intensity measurement, 
following |Wardle & Kronber^ ( |1974[). Detail s about the polariza¬ 
tion calibration can be found in|Surcis et al.|(|20T3]l. 


3. Results 


In the HMC center of G023.01-00.41, the milli-arcsecond distri¬ 
bution of the 6.7 GHz CH 3 OH maser cloudlets has not changed 
over the five years spanned by our observations. Although at 
some maser velocities the overall flux density has smoothly 
changed in time (Eig.[T]), this variation affects the flux density 
of individual cloudlets, whereas the overall maser distribution is 


preserved. Among the eighty maser cloudlets detected by Sanna 
|et al.| ( | 2010 b| l, only about one third shows linearly polarized 
emission, and their properties are listed in Tab.[^ Given that 
we know only the orientation of the linear polarization vectors 


* The European VLBI Network is a joint facility of European, Chinese, 
South African and other radio astronomy institutes funded by their na¬ 
tional research councils. 

^ http://www.aoc.nrao.edu/ smyers/evlapolcal/polcal_master.html 
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Fig. 2. Comparison of the velocity and polarization vectors across the 
CH 3 OH maser distribution in G023.01-00.41. We show the distribu¬ 
tion of the 6.7 GHz maser cloudlets (dots) within a field-of-view of 
about 2000 AU, following [Sanna et al.| l |2010b| their Figure 6 ). Maser 
velocities along the line-of-sight are color coded according to the right- 
hand Vlsr scale. Lines drawn on individual maser cloudlets indicate the 
orientation of the velocity (red) and polarization (black) fields locally, 
as obtained from VLBI measurements. For clarity, velocity vectors are 
only shown for cloudlets with an associated polarization vector. For iso¬ 
lated maser cloudlets, numbers correspond to the maser labels given in 
Tab.[^ The star marks the peak position of the high-excitation, CH 3 OH, 
thermal line detected by jSanna et al.|p014^ , which is the origin of the 
reference system (e.g., upper right corner) used in Fig.[^and|^ 


(UpoOj while their direction is undefined, we folded these val¬ 
ues in the range -90°<;ifpoi<90°. The linear polarization fraction 
detected among the maser cloudlets ranges between 0 . 6 % and 
9.2%, with an average value of 2%. The orientation of the linear 
polarization vectors is superposed on the CH 3 OH maser distri¬ 
bution in Eig.[^ For cloudlets with detected linear polarization, 
we also draw the direction (Uvei) of the velocity vectors for com¬ 
parison. 

In Fig.[^ we study the distribution of the polarization vec¬ 
tors orientation with respect to the position of the CH 3 OH maser 
cloudlets. In Fig.j^, we plot the minimum difference (Upoi-Avei) 
between the orientation of the polarization and velocity vectors 
as a function of the sky position angle of each cloudlet. This po¬ 
sition angle is measured east of north with respect to the HMC 
center (star symbol in Fig.|^, which is defined as the peak posi¬ 
tion of the high-excitation, CH 3 OH, thermal emission at the sys¬ 
temic velocity of the core ( jSanna et al.|2014[ ). This position rep¬ 
resents the current best estimate of the YSO position. The Xpoi- 
distribution has a weighted average of 58°, and a weighted 
dispersion of +18° (grey area in Fig[^). In Fig.[^ andc, maser 
positions are projected along the two orthogonal axes defined by 
the direction of the molecular outflow (+58°), and that of the 
elongated HMC (-32°), with the origin at the HMC center. As¬ 
suming a perfect symmetry of the gas dynamics with respect to 
these two axes, we produce a mirror image of the linear polariza¬ 
tion vectors of each cloudlet on a single quadrant (Xpo\)- For the 
x'p^^ values, we also try to solve the ambiguity of +180°, assum- 
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ing that cloudlets close in space would show a smooth change 
of the polarization (and magnetic) field with position. A poste¬ 
riori, this criterion is found to minimize the difference between 
the direction of the linear polarization vectors and that of the 
corresponding velocity vectors. In Fig.[^ andc, we identify two 
regions, labeled Reg. 1 and 2, where changes smoothly with 
the projected distance of the maser cloudlets. These regions are 
fitted by a linear slope of about 0.2° AU ', and correspond to 
nearby cloudlets which also show a smooth variation of the ve¬ 
locity field (Fig[^. A third region labeled Reg. 3, shows a vari¬ 
ation of;ifp^j by more than 100 ° over a small range of projected 
distances. This region corresponds to cloudlets with the shortest 
projected distances, and will be discussed further in Sect.j^ 


4. Discussion 


To infer the local magnetic field orientation, we ran the radiative 
transfer model by Vlemmings et al. ( 2010[ ) for each cloudlet with 
detected linearly polarized emission (last columns of Tab.[2l. Ac¬ 
cording to the output parameters of this modeling (0» 55°), the 
magnetic field orientation is perpendicular to the polarization 
vectors for all maser components. This information is used in 
Fig. to plot the local, magnetic field orientation by rotating the 

values by 90°. More details about the radiative transfer mod¬ 
eling used, can be found in |Surcis et al.| ( |20T3) l . 

In Fig .|4] we make use of the outflow/disk geometry de¬ 
scribed in Sanna et al. p014| l, to give a complete picture of the 
velocity and magnetic fields within 1000 AU from the FIMC 
center. In this plot, we produce a mirror image of all the mea¬ 
surements of velocity and magnetic field vectors obtained from 
the 6.7 GHz, CH 3 OH, maser cloudlets, as if they were sampling 
a single quadrant defined by the outflow direction and the disk 
plane. This picture holds under the assumption that the gas dy¬ 
namics shows a symmetric behavior with respect to the outflow 
axis and the disk plane. Given the uncertainty of 30° on the 
disk inclination ( |Sanna et al. 2014| l, one should keep in mind 
that the maser cloudlets might be closer to the disk plane than 
they appear. The high-density molecular tracers observed toward 
G023.01-00.41 show a fairly constant ratio of ~2 between the 
major and minor axis of the HMC. In Fig.|^ this ratio is inter¬ 
preted as if it was due to a flared disk with a semi opening angle 
of 30° (dark grey area). We also mark a central region along the 
outflow axis which is devoid of CH 3 OH maser emission (light 
grey area). 

In an attempt to derive a continuous, magnetic field morphol¬ 
ogy, which reproduces the local maser measurements, we con¬ 
sidered those regions (Reg. 1 and 2) showing a smooth change 
of the polarization vectors (and magnetic field) orientation with 
the maser cloudlets position. At a first order, these slopes have 
been approximated by a linear fit as shown in Fig.|3|(;t'p^j=/(.r)). 
We can then integrate the tangent of [f(x) + 90q, in order to 
derive the families of curves which best fit the local, magnetic 
field orientation at the maser cloudlets position. These curves 
give a first order representation of the local morphology of the 
magnetic field lines (black dotted lines in Fig|^. 

The velocity field traced by maser cloudlets belonging to 
Reg. 1 and 2, provides a consistent picture of gas outflowing 
from the HMC center along the magnetic field lines. In Reg. 1, 
for small heights over the disk plane (<400 AU), the velocity 
field well accommodates with the magnetic field lines, start¬ 
ing from projected distances of 600 AU up to about llOOAU. 
Further away, the velocity field mainly expands parallel to the 
disk plane, and shows the highest maser velocities (~ 10 km s '). 
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Fig. 3. Analysis of the polarization vectors orientation (Xpoi) as a func¬ 
tion of the maser cloudlet position (see Sect.|^. a) plot of the difference 
(and uncertainty) between the orientations of the velocity ixva) and po¬ 
larization (Xpo]) vectors as a function of the sky position angle (P.A.) 
of each cloudlet (defined in Sect.[^. Cloudlets with P.A. in the range 
90°-270° have been folded in the range -90°<P.A.<90°. The grey area 
marks the boundary of the weighted standard deviation (± 1 cr) for the 
Apoi-Avei values, b) plot of x'j,oi as a function of the projected distance 
along the equatorial plane of the HMC (-32°). x'^^x gives the values of 
Xpoi as measured counterclockwise starting from a sky position angle 
of -32°. In this panel, only points with offsets less than 400 AU from 
the equatorial plane have been plotted. For each value, the polar¬ 
ization fraction is quantified according to the wedge on the right-hand 
side. The uncertainty of the x'^^i measurements is of the same order as 
the marker size. The different regions identified in Sect.j^are also indi¬ 
cated. c) similar to panel b) with projected distances along the outflow 
axis (58°). Data points plotted in panel b) are not drawn. 


On the other hand, maser cloudlets belonging to Reg. 2 and up¬ 
wards, expand and get collimated in the direction of the out¬ 
flow axis. In particular, as one proceeds upwards along the mag¬ 
netic field lines, and closer to the outflow axis, both the veloc¬ 
ity and magnetic field vectors independently undergo a turn of 
90°. This feature may be interpreted as a result of the complex 
gas dynamics where both a slow and fast velocity component 
exist (e.g., Seifried et al.||201^ their Fig. 12). Indeed, at about 
2000 AU along the outflow axis (not shown in Fig.[^, a shock 
front of dense gas traced by H 2 O masers shows gas velocities of 
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Fig. 4. Combined plot of the velocity and mag¬ 
netic field measurements obtained from the 6.7 GHz, 
CH3OH, maser cloudlets (dots) in G023.01-00.41. 
The axes are aligned with the outflow and disk direc¬ 
tions, as defined by the sub-arcsecond SMA obser¬ 
vations in the 1 mm band from |Sanna et al.|P014| l. 
The origin is set at the HMC center, as defined in 
Sect.|^ (star symbol in Fig.[^. The dark grey area 
marks the region where emission from a flared disk 
with axes ratio of 1 /2 is expected. The light grey area 
marks the region close to the outflow axis which is 
devoid of maser emission. For each maser cloudlet, 
the local directions of the velocity vector (red solid 
arrows) and/or that of the magnetic field vector (black 
solid lines) are reported where available. Empty dots 
are used for cloudlets which do not show polarized 
emission. The local average of both the magnitude 
and direction of the velocity field is indicated by red 
dotted arrows. These values were derived by group¬ 
ing close by features with similar velocity directions. 
The different regions identified through the polariza¬ 
tion vectors analysis (Fig.|^, and labeled from 1 to 
3, are indicated. For regions showing a smooth vari¬ 
ation of the polarization (and magnetic) field vectors 
with the maser cloudlets position, we plot the fami¬ 
lies of curves (black dotted lines) which best match 
the tangent to the local, magnetic field measurements 
(see Sect.|^. The pattern of these families of curves 
approximates the spatial morphology of the magnetic 
field lines locally. The direction of the velocity vec¬ 
tors is on average tilted by 30° with respect to the 
local, magnetic field direction. 


20km s ^( Sanna et al.||2010b their Fig. 5b), three times higher 
than those traced by the CH 3 OH gas in Reg. 2. 


Methanol maser cloudlets belonging to Reg. 3 have projected 
distances of less than about 500 AU (dashed box in Fig. 1 ^, both 
along the disk plane and the outflow axis. In this inner region, the 
velocity field is composed of (at least) two different motions, 1 ) 
an inflowing motion closer to the disk plane (y-offset < 200 AU), 
and 2) an upward motion for higher offsets. The magnetic field 
pattern of this region is more complex, and we did not attempt to 
reproduce the magnetic field lines. Still, nearby cloudlets show a 
similar orientation of the magnetic field vectors, and confirm the 
accuracy of our measurements. 


Interestingly, Reg. 3 corresponds to a diffuse halo emis¬ 
sion from strong CH 3 OH masers ( |Sanna et al.||2010b see their 
Fig. 6 b), which are likely saturated due to the vicinity of the cen¬ 
tral IR source. This ridge of extended emission is significantly 
elongated in the direction of the outflow axis, in agreement with 
the average orientation of the magnetic field vectors between y- 
offset~300 and 400 AU, and the upward motions detected there. 
This evidence makes us speculate that Reg. 3 may trace the outer 
launching region of the primary outflow, in agreement with re¬ 
cent MHD simulations by |Seifried et al.| ( |2012| their Fig. 5). 
Furthermore, we make use of the inward stream of gas close 
to the disk plane, to obtain an estimate of the mass inflow rate, 
Min = (5.0 X 10“^Moyr"')RjgpVion8. In this formula, Rioo, 
vio, and ng are the mean radius of the inward stream in units of 
lOOAU, its velocity in units of 10kms ', and the volume den¬ 
sity of molecular hydrogen in units of lO^cm^^, respectively. 
We take into account that the inward stream of gas is confined 
within an angle of 60° either side of the YSO, and allow for gas 
densities as high as 10*'cm“^, above which the Class ii CH 3 OH 
masers start to be quenched ( Cragg et al.|2005| ). Noticeably, we 
find no maser detection at closer distances to the HMC center. 
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At an average distance of 300 AU, the inward stream of gas, 
flowing at a velocity of 5kms“*, brings a mass inflow rate of 
2x 10-4 Mo yr-'. 


We finally consider the average misalignment of 60° ob¬ 
served between the velocity and polarization vectors (Fig.[^), 
which translates to 30° between the velocity and magnetic field 
vectors (Fig.|^. If gas and magnetic field were fully coupled, we 
would observe gas flowing along the magnetic field lines, and 
expect an average misalignment close to zero. Given that the ob¬ 
served misalignment appears randomly distributed (1 cr = +18°) 
about 30°, we model this effect as it was due to a random (turbu¬ 
lent) velocity component, which adds to ordered velocity vectors 
aligned with the magnetic field lines. By considering a velocity 
vector of 7 km s^*, as averaged across the whole region, we esti¬ 
mate this turbulent contribution to be of the order of 3.5 km s '. 
This value is very similar to the velocity dispersion (4-5 km s^') 
derived from the CH 3 CN linewidth in the inner 3000 AU from 
the HMC center ( Sanna et al.||2014) l, which fairly supports our 
estimate. 
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Fig. 1. EfFelsberg total-power spectra toward G023.01-00.41, obtained from EVN observations at C-band, over 5 yr (from Sanna et al.|2010b and 
current measurements). Observing dates are indicated on the top right of the plot. The dotted vertical line marks the systemic velocity (Vj^j) of the 
HMC, as inferred from CH3CN measurements. 
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Table 1. Parameters of 6.7 GHz methanol maser cloudlets with detected linear polarization. 


Feature 

Vlsr 

Fpeak 

Pr 

/^pol 

Xy&I 

Radiative Transfer Model 
AVi TbAQ 0 

# 

(kms-') 

(Jybeam”*) 

(%) 

(°) 

n 

(kms"') (logKsr) (°) 


Northern Region - Dec. (J2000) > -9:00:38.30 


1 

74.76 

105.50 

4.0 

± 

0.1 

-61 

± 

2 

+22 + 14 

1.5!°'2 

2 

75.55 

34.52 

1.1 

± 

0.1 

-14 

± 

1 


2 6 ^ ■ 

3 

74.32 

44.47 

2.2 

± 

0.3 

-22 

± 

4 

+ 108+19 

1 4+0-2 
^ -^-0.4 

4 

74.54 

97.84 

1.1 

± 

0.1 

+51 

± 

1 

+5 + 7 

9 9 + 0.2 
-0 3 

5 

74.81 

45.39 

0.9 

± 

0.2 

+81 

± 

2 



7 

75.86 

38.56 

1.1 

± 

0.1 

+76 

± 

1 


2-3+! 

11 

72.70 

30.90 

2.2 

± 

0.1 

+36 

± 

5 

+ 140+14 

2 

^■^-04 

13 

74.41 

8.24 

4.6 

± 

0.8 

+36 

± 

2 


1 l +°'2 

14 

74.32 

19.42 

3.3 

± 

0.3 

+27 

± 

1 


1 4 + 0:2 

■ -0 5 

15 

74.67 

11.37 

0.9 

± 

0.3 

0 

± 

9 

+96 + 8 

1 4 + 0.2 
^ ■ -0 2 

18 

73.75 

9.48 

1.7 

± 

0.3 

+ 17 

± 

5 

+75 + 7 

l-4!nl 

19 

73.49 

10.69 

1.4 

± 

0.1 

-17 

± 

1 

-160+10 

2 . 2^1 

21 

76.60 

10.09 

1.4 

± 

0.3 

-10 

± 

3 

+97 + 6 

2 4^ 
- 0.3 

24 

73.93 

4.03 

1.5 

± 

0.3 

Oi 

12 


< 0.5 

26 

73.18 

11.29 

0.6 

± 

0.2 

-88 

± 

6 

-24 + 10 

1 3^°' 

31 

73.27 

4.32 

3.8 

± 

0.4 

+33 

± 

2 

+ 147 + 38 

1 6 +°-^ 
-0 5 

65 

78.71 

0.82 
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Notes. Labels in Col. 1 are used to associate the current measurements with cloudlets identified in Tabled of |Sanna et al.| ( |2bl0b^ . Columns2 
and 3 report the LSR velocity and brightness of the brightest spot of each cloudlet on 2011 May 29. Columns 4 and 5 give the measured linear 
polarization fraction and the position angle (east of north) of the linear polarization vectors for each cloudlet, respectively. Column 6 gives the 
position angle of the velocity vectors as measured from |Sanna et al.|j2010b| l, for comparison with column 5. For each maser, Cols. 7, 8, and 9 give 
the model results for the emerging brightness temperature, the intrinsic thermal linewidth, and the angle between the magnetic field and the maser 
propagation direction, re^ectively. 

New feature, see Fig.^ 













